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Wc find two systems via holography that exhibit quantum Berezinskii-Kosterlitz-Thouless (BKT) 
phase transitions. The first is the ABJM theory with flavor and the second is a flavored (1,1) little 
string theory. In each case the transition occurs at nonzero density and magnetic field. The BKT 
transition in the little string theory is the first example of a quantum BKT transition in (3+1) 
dimensions. As in the "original" holographic BKT transition in the D3/D5 system, the exponential 
scaling is destroyed at any nonzero temperature and the transition becomes second order. Along the 
way we construct holographic renormalization for probe branes in the ABJM theory and propose 
a scheme for the little string theory. Finally, we obtain the embeddings and (half of) the meson 
spectrum in the ABJM theory with massive flavor. 



I. INTRODUCTION 

Holography [TrEj has become a cornerstone in the 
study of strongly interacting theories. Its use has ex- 
tended far beyond its strict range of applicability, afford- 
ing simple geometric pictures for the physics of confine- 
ment [4], chiral symmetry breaking [5], and phase tran- 
sitions [6 . Indeed, holography has found many more 
applications to heavy- ion physics, including the study of 
transport in strongly coupled plasmas [7] and the energy 
loss of hard partons [5J |S] . More recently there has been 
a flood of work realizing condensed matter phenomena 
via holography. These have included the Fermi gas at 
unitarity [101 [TT], the quantum Hall effect P2HH], su- 
pcrfluids [T5J [TB] , non-Fermi liquids [T7H2"U] , and quan- 
tum critical points [2"TM24| . Reviews of many of these 
applications can be found in [2"S"H2"7] . 

The study of critical phenomena is of central impor- 
tance in the condensed matter community. In addition 
to its many storied successes [55] there are some unsolved 
problems surrounding the role of strongly interacting 
quantum critical points and strange metals [29] . These 
are relevant in candidates for the theory of high T c su- 
perconductors. It is a driving hope that holography can 
successfully describe these systems. A first step toward 
this goal is to identify the classes of phase transitions that 
are natural in holographic theories. Unfortunately, most 
of these transitions are rather boring from this point of 
view: they are usually first-order 4, 6j 30]|3T] or second- 
order with mean- field exponents 32 34J . This mean- field 
scaling is actually expected rather than surprising; there 
is a large N parameter that suppresses quantum fluctu- 
ations in both the field [35] and gravitational theories. 
This represents a challenge to the study of more interest- 
ing phase transitions at large N: in additional to identi- 
fying them [3H1 [37] , we should also understand why they 
exist in the first place in the sense of [38] . 

The first holographic example of a continuous non- 
second-order transition embedded in string theory was 



* kristanj9u.washington.edu 



found in the D3/D5 system [22:. There is a novel chi- 
ral quantum transition in that system with the scaling of 
the Bcrczinskii-Kosterlitz-Thouless (BKT) phase transi- 
tion [351 EH]- Recall that transitions of the BKT type 
are between disordered and Quasi-ordered phases in two 
dimensions. BKT transitions are rather special: their 
existence is intimately related to the Coleman-Mermin- 
Wagner theorem [41"] - |4"3] . The two-point function of the 
order parameter in the disordered phase exhibits a cor- 
relation length that scales with exp(c/y/T — T c ) near the 
critical temperature T c [441 . Similarly, the free energy of 
the quasi-ordered phase differs from that of the ordered 
by an amount that scales as exp(— cj \JT C — T). 

The transition in the D3/D5 system is novel in that 
it has BKT scaling in an ordered phase. For this reason 
we termed it a holographic BKT transition in [52], since 
it happens in a different context from the original BKT 
transition. Before going on, we note two more interesting 
features of the critical point in the D3/D5 system. First, 
the exponential scaling of the order parameter is lost for 
any nonzero temperature, for which the chiral transition 
becomes second-order with mean-field exponents [22]l45]. 
Second, there are relics of the BKT scaling at nonzero 
temperature: the critical temperature of the transition 
scales exponentially [2"3] . 

We do not often stumble across new types of phase 
transitions. As a result, we should ask ourselves a num- 
ber of questions: (i.) Can we identify other field theories 
that realize holographic BKT transitions, particularly in 
theories without a gravitiational dual? (ii.) What is the 
field theory mechanism that triggers these transitions? 
(iii.) Can we identify an effective theory near critical- 
ity? In this work we attack the first question by identi- 
fying more examples. However, our goal is not to simply 
enumerate a list of theories with holographic BKT transi- 
tions. Rather, we seek to develop a classification program 
toward the end of answering questions (ii.) and (iii.). 

In this work we study flavored ABJM theory 46] and a 
flavored (1, 1) little string theory [17] in detail. These are 
supersymmetric theories that we study holographically. 
In particular, wc identify a holographic BKT transition 
in both systems at nonzero density and magnetic field. 
The gravitational description of these systems is given in 
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FIG. 1. A plot of the condensate in the flavored little string 
theory as a function of magnetic field at zero and finite tem- 
perature near the zero-temperature transition. The dashed 
black line indicates zero-temperatur e num erical data and the 
solid blue line our prediction, Eq. (3.351. The color dashed 



curves represent numerical data at background entropy densi- 
ties of s = 10~ 24 /(2tt) 4 ^J? (left) and 10" 22 /(2vr) 4 ^ii (right). 
At any nonzero temperature, the condensate scales with a 
mean-field exponent near the transition and then asymptotes 
to the BKT scaling at large magnetic field. 



terms of a probe brane minimizing its worldvolume action 
in a fixed geometry. In Refs. [ST] [25] we studied similar 
probes describing (3 + 1) and (2 + 1) dimensional field 
theories. The (3 + 1) dimensional system was also stud- 
ied in |48j . There the competition between a nonzero 
density and magnetic field gave rise to chiral quantum 
transitions: the transition in the (3 + 1) dimensional the- 
ory was second-order with mean-field scaling, while the 
same system in (2 + 1) dimensions exhibits exponential 
scaling. In this work, the flavored ABJM theory also lives 
in (2 + 1) dimensions but the flavored little string the- 
ory is the first example of a theory in (3 + 1) dimensions 
with BKT scaling. For dramatic effect and summary, we 
have included a plot of the order parameter in the broken 
phase of that theory in Fig. [T] 

The gravitational mechanism for the holographic BKT 
transition is clear. It occurs due to the violation of the 
Breitenlohner-Freedman (BF) bound [49) in the infrared 
region by the scalar field dual to an order parameter. On 
general grounds presented in [50] this violation was ex- 
pected to produce BKT scaling. However, there are only 
a handful of known settings where the BF bound is vio- 
lated in a controlled setting. The first was the D3/D5 sys- 
tem, the second in extremal asymptotically AdS4 dyonic 
black holes [33], and the third in this wor 10 While this 



mechanism seems similar to that of a second-order Lan- 
dau phase transition, the bulk field actually represents 
an infinite tower of field theory modes. At the transi- 
tion, an infinite number of these modes destabilize. This 
situation is quite unnatural within the Landau theory. 
To us these results suggest that we need a new paradigm 
of phase transitions to describe the infrared physics of 
these systems. 

While we have not yet found such a paradigm, we have 
identified a picture for some sufficient conditions to trig- 
ger a holographic BKT transition. We do this by consid- 
ering a large class of theories without holographic BKT 
transitions. In particular, we add flavor to the field the- 
ories dual to the near-horizon p-brane geometries. These 
are interacting theories in (p + 1) dimensions with run- 
ning couplings (except for p — 3, corresponding to Af = 4 
super Yang-Mills (SYM) theory) but with a generalized 
conformal symmetry. Operators typically carry different 
dimensions than their canonical ones, a feat accomplished 
by redefining them with powers of the couplings. After 
turning on a density and magnetic field, we find that the 
only theories that realize holographic BKT transitions 
are the D3/D5 system and the flavored little string theory 
of this work. Combined with flavored ABJM, these are 
also the only theories for which the generalized dimension 
of the density and magnetic field equal each other. We do 
not think that this is a coincidence. From these results 
we conjecture that a holographic BKT transition can be 
triggered in a generalized conformal theory after adding 
equal-dimension control parameters with some other con- 
straints. See our discussion in Sec. HVI for more details. 

There are a few other cute results that we establish 
while studying these BKT transitions. The first is in the 
ABJM theory. In order to properly study chiral sym- 
metry breaking in that theory, we find it convenient to 
obtain the brane embeddings dual to the theory with 
supersymmetric massive flavor. While we do not explic- 
itly verify kappa-symmetry, we perform a few consistency 
checks that these branes are indeed supersymmetric. Af- 
ter holographically renormalizing the bulk theory, we find 
that the free energy and condensates of the field theory 
vanish. For another check, we obtain half of the meson 
masses of the dual theory. These are suppressed from the 
quark mass by a factor of the 't Hooft coupling of the 
ABJM theory, m mcBon ~ "iquark/vA. We also find hints 
of a broken 5*0(5) invariance in the meson masses, much 
like in the D3/D7 system [55]. Finally, in the course of 
regularizing the little string theory we also holographi- 
cally rcnormalize the D3/D3 system at nonzero density. 
The most interesting result here is a new Weyl anomaly 
of the dual theory proportional to /i 2 , where /i is the 
chemical potential for baryon density. 

The outline of this work follows. In Sec. [IT] we describe 
the probe brane setups we use, beginning with flavored 



1 While we were finishing this work, S. Pal released a paper |51l 
that also studies holographic BKT transitions in the Dp/Dq sys- 
tems. Their results at nonzero density and magnetic field agree 
with ours. 



ABJM theory in Sec. |II A| First we review the ABJM the- 
ory and the probe branes dual to the addition of massless 
flavor. We go on to construct the embeddings dual to su- 
persymmetric massive flavor and compute some the me- 
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son masses of that theory. Next, we review the (1,1) little 
string theory and the addition of flavor in Sec. |II B| In 
Sec. |II C| we holographically renormalize probe branes in 
the AB JM theory and propose a scheme for a holographic 
map in the little string theory. The bulk of our results are 
in Sec. |III[ where we analytically and numerically identify 
the holographic BKT transitions in these theories. We 
also perform numerics for the little string theory, mea- 
suring the fate of exponential scaling at nonzero temper- 
ature in Sec. |HIB~3"l Finally, we discuss our results in the 
context of the Dp/Dg systems and generalized conformal 
symmetry in Sec. |IV| 



II. THE BRANE SETUPS 

A. ABJM with flavor 

The geometry supported by a large stack of N M2 
branes in M-theory probing the orbifold singularity of 
C 4 /Zfc has a near-horizon limit of AdS4 x S 7 /Zfc, 
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(R/hi) 6 = 32ir 2 kN, 



(2.1) 
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with kN units of 4-form flux through the AdS4 
tor. The field theory living on the branes is Af = 6 
U(N) x U(N) superconformal Chern-Simons theory at 
level k gB]. This is the celebrated ABJM theory. At 
levels k = 1,2, the theory has enhanced Af = 8 super- 
symmetry [46J [53] . Moreover, the theory at N = 2 is 
identical to the Bagger-Lambert theory [Ml |SS] proposed 
a few years ago as the worldvolume theory on a stack of 
M2 branes. 

The orbifolded circle inside § 7 /Zfc has a size R/k which 
is small when N <C fc 5 . In this limit the M2 branes are 
well-described by the compactification of AdS4 x § 7 /Zfc 
on this circle, i.e. by type IIA string theory on AdS4 x 
CP 3 . In string frame, this background is 



o 20 _ 



R 3 
4A- 

k 3 ' 



(5AdS 4 + ^9CF 3 )i 

F 4 = |i? 3 vol(AdS 4 ), F 2 = kJ, (2.2) 



where we have set a' = 1. The two-form F2 appears in 
the lOd description because the reduction is on a non- 
trivially fibered circle inside § 7 /Zj,; the connection asso- 
ciated with the fibration is the potential for the Kahler 
form on CP 3 , J. This background was first identified 
in [56] as the dimensional reduction of the AdS4 x S 7 
solution of lid supergravity on the diagonal circle in S 7 . 

We therefore have a IIA description of the ABJM the- 
ory at large k and N when k 5 ^> N and when the lOd 
radius of curvature is large, i.e. R 3 /4k oc yWfk 3> 1. 
Most of the applications of this duality have centered 
on the conjectured integrability of the ABJM theory by 
studying the dimensions of large spin operators 



Additionally, there has been some work done adding fla- 
vors to the ABJM theory [65H68] . This addition is re- 
alized in the bulk by adding a small number of probe 
D-branes [55] that wrap some cycles inside the CP 3 . 

It is a short exercise to show that adding a single 
fundamental hypermultiplet to one of the gauge groups 
breaks the R-symmetry of the theory to 50(3) [6ST - 
157] and so the supersymmetry to Af = 3. There is 
an additional SO(3) x global symmetry preserved when 
the flavor is massless so that the global symmetry is 
SO(3) R x S0(3) x x U(1) B , where the last U(l) sym- 
metry is a baryon number symmetry. By supersymme- 
try the dual brane setup should have a single D6 brane 
wrapping a 3-cycle inside of CP 3 . Additionally, this cycle 
should have a Z 2 fundamental group, so that there can be 
two different Wilson lines on the brane corresponding to 
adding the flavor to one gauge group or the other [631170] , 

Up to SU(4) rotations of CP 3 there is a single SO (4) 
invariant 3-cycle, EP 3 . This cycle also has a Z2 fun- 
damental group, making it a consistent candidate for 
the correct flavor brane embedding. Indeed, a D6 brane 
wrapping AdS4 x EP 3 preserves twelve supercharges of 
supersymmetry [65 . In the probe limit massless flavored 
ABJM theory is a superconformal Af = 3 theory, so this 
is the correct dual. 

However, our ultimate goal is to study SO(3) x chiral- 
symmctry breaking in this theory. The embeddings dual 
to the chirally broken phase of this theory will deviate 
from the supersymmetric embedding. The deformation 
should be in the same direction as that of the embeddings 
dual to massive flavor, which have not yet been obtained. 
Thus, in order to correctly study SO(3) x breaking we 
elect to study the addition of massive flavor. This is 
actually not that hard. 

In the bulk we will have a dual flavor brane with a non- 
trivial embedding and zero field strength for the U(l) 
gauge field on the brane. The non-trivial part of this 
problem is that these embeddings will deform the EP 3 as 
a function of the radial coordinate. We will therefore look 
more carefully at the embedding of EP 3 inside of CP 3 . 
Following [55], we find it instructive to construct CP 3 
inside of C 4 . We will use this construction to describe 
EP 3 and its deformations. 



1. Massive flavor 

We begin with N M2 branes on C 4 , which we 
parametrize with complex coordinates Zi , i — 1 — 4 as 

9 1 - xi+4>i , 9 1 - xi-<Pi 

z\ = r cos £ cos — e 2 , Z2 = r cos f sin— e 2 , 

. 9 2 . . 62 3 *2--»2 

Z3 = r sin£ cos — e 2 , 2:4 = r sin£sin — e ^.3) 

where r is the radial coordinate of C 4 and the angular 
variables have the domains £ € [0,5],^ € [0, 7r],Xi € 
[0, 4w),(pi £ [0, 27t). Defining new angular coordinates y 
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and ip by 



angular part S 7 /Zfe as 



Xi = 2y + ip, X 2 = 2y- ip, 



(2.4) 



we see that the diagonal circle of C 4 is given by y £ [0, 2ir) 
along which the orbifold acts. The orbifold simply 
identifies y ~ y + ^ , so that we write the metric of the 
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which simply says that S 7 /Zfc can be written as a Hopf 
fibration over CP 3 with connection A. In these coordi- 
nates we have 

2A = cos 2£cfa/> + cos 2 £ cos 9id(pi+sm 2 {; cos Qid<pi, (2.6) 

and 
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The connection A is the potential for the Kahler form on 
CP 3 , dA = J. 

The MP 3 inside of CP 3 we consider is invariant under 
the simultaneous SU(2) rotations of (21,22) and (23,24) 
(the bulk realization of the 5*0(3)r, symmetry) and is 
given by 
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The induced metric on this space is 
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where the coordinates have domains 6 € [0, 7r],0 € 
[0, 2-71"), ?/> G [0, 2ir). This is indeed the metric for a round 
unit RP 3 . It is also invariant under a second 577(2) 
that simultaneously rotates (21,24) and (z 3 ,2 2 ). This 
5?7(2)2 ~ 5*0(3) x is the bulk realization of the extra 
chiral symmetry of the massless theory. 

Of course this is not the only 3-cycle invariant under 
the SO(3)b,- There is an obvious SO(3)r invariant 3- 
cycle 
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which becomes RP for £0 = tt/4. The ^-direction is 
transverse to this cycle as are two other directions re- 
lated by S*0(3) x transformations. The fluctuations of the 
cycle under these three directions transform as a triplet 
under SO(3) x . The dual operator is the SO(3) x triplet 
hypermultiplet mass operator ipip + SUSY [65] . We find 
it convenient to consider brane embeddings where the 
brane slips off of MP 3 in the ^-direction. 

There is a useful change-of-coordinates for this prob- 
lem. Define y, p by 



7r 1 y 
t = — I — arctan — , 
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In these coordinates the Poincare horizon of AdS4 is lo- 
cated at y = p = and the boundary of AdS4 at either 



y, p — > 00. Then a D6 brane that wraps C3 with £ = £(r) 
can be described with an embedding y = y(p) and an 
induced metric 
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where the internal space has a metric 
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We recognize this space as a squashed unit 
squashing parameter y/p, 



9c 



\ y - 2 (do 2 

4 p A 



sin 



(2.13) 
with a 
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When y = the internal space is a round RP 3 and as 
y/p — > 00 the brane ends smoothly as detP[g] — > 0. 
This change-of-coordinates is similar to the one in [22] 
that has made many computations in probe brane sys- 
tems tractable. Notably, this squashing only preserves 
the U(l) in SO(3) x that rotates ip. The global symme- 
try of the dual theory is then 5*0(3) R x U(l) x x U(1) B . 

Before obtaining the embeddings dual to massive fla- 
vor, we need one more ingredient. There is a non-trivial 
Chern-Simons term on the D6 branes that must be in- 
cluded in order to obtain the correct equations of motion 
for the embedding, the integral of C7 on the worldvolumc. 
The pullback of the CV to the brane is [HS] 

R 9 

P[C 7 ] = Tr^r 2 sin 9dx° A dx 1 A dx 2 AdrAdipAd0A dep, 

(2.15) 



and the brane action is 



T 6 / P[C 7 ] (2.16) 



in the absence of any worldvolume field strength. This 
Lagrangian for the embedding function y is rather com- 
plicated. Integrating over the internal space and dividing 
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by the volume of the field theory directions , it is pro- 
portional to 

£ 6 oc -pTp^yTTF + f {P r^\ ■ (2.17) 

The equation of motion for this action has a remark- 
ably simple solution, y =constant. We claim that these 
embeddings describe the ABJM theory with massive su- 
persymmetric flavor: we have consistent supergravity so- 
lutions dual to the field theory where we have turned on 
a source for the hypermultiplct mass operator. Indeed, 
as we show in Sec. 
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(i.) there is a renormalization 
scheme under which the dual extremum has zero free en- 
ergy and (ii.) the expectation value of the field dual to y 
vanishes. Both of these results are important consistency 
checks with the supersymmetry of the dual theory. Yet 
another check, the meson spectrum, is mostly computed 
in Appendix [A] and is reviewed below. 



2. The meson spectrum 

We obtain the meson spectrum of flavored ABJM the- 
ory by studying the fluctuations of worldvolume fields 
around the embedding y — m. In order to be a small 
fluctuation these fluctuations must be small and well- 
behaved everywhere on the brane, particularly near p = 
0. This is a boundary condition on the fluctuations. 
Fourier-transforming in the field theory directions x^, the 
fluctuations will only obey a normalization condition near 
the AdS4 boundary for particular momenta —k 2 . These 
are the masses (-squared) of the mesons of the dual field 
theory. 

For our choice of parametrization there are three 
classes of worldvolume fields, (i.) the transvese scalar 
y and those related by SO(3) x symmetry, (ii.) the U(l) 
gauge field A, and (iii.) a neutral fermion 'J. As we show 
in Appendix |Aj the fluctuations of the gauge field can be 
broken up further into three types, a vector and two dif- 
ferent types of scalars. We elect to study the bosonic 
mesons on the presumption that supersymmetry relates 
their masses to the fermionic mesons. 

The supercharges of the theory fill a triplet of the 
50(3)r symmetry. These commute with the S*0(3) x chi- 
ral symmetry, which is broken to U(l) x by turning on a 
mass. All states in a given supermultiplet then have the 
same U(l) x charge. Our mesonic multiplets can then 
be obtained by acting with the supercharges on a scalar 
state with spin I under the SO(3)r, that is annihilated 



by the Q's. Before going on, we will use the notation 
(2' if' ir) to denote the quantum numbers of a meson, 
where are S*0(3)r and SO(3) x spins respectively and 
?Y is the U(l) x charge. Also, j, j,' , and m' must be even, 
which is related in the bulk to the fact that the internal 
space C3 is a Z2 identification of a squashed 3-sphere. 

We are not able to obtain the full meson spectrum of 
the theory. We are, however, able to obtain several sec- 



tors. We describe their computation in Appendix [A"| and 
simply summarize here. First, we solve the fluctuation 
spectrum of the transverse scalar y and find the meson 
masses 



fc 2 = — (2+j + \m'\+4n)(i+j + \m'\+4n), (2.18) 



in the |, ,J ^j represenation for any positive integer 

n. The operators dual to these fields have dimension 
2 + |. The j — operator is the hypermultiplet mass 
operator ^/> + SUSY, which has dimension 2 as expected. 
The other transverse scalars are related to y by SO(3) x 
symmetry but they have U(l) x charges ±1. The masses 
of these dual mesons are then given by Eq. ( 2.18[ ), but 
with m! shifted by =f2. The net effect is that for |m'| < j 
there are three sets of degenerate meson spectra with 
masses Eq. (2.181 and that there are extra mesons with 
U(l) x charge ± (| + l) and masses 



k 2 = m 2 (l + j + 2n)(2 + j + 2n). 



(2.19) 



We are also able to obtain the masses of the vector 
mesons, dual to eigenmodes of a type of gauge field. 
These obey the same equation as the transverse scalar 
and so have the same spectrum as Eq. (2.181. These also 

live in the |, representation and are related to 

operators of dimensions A = 2 + | . 

There are a few sectors that we were not able to obtain. 
The first were another set of scalars in the f|; |, 
representation dual to another fluctuation of the gauge 
field. We did, however, study this sector numerically and 
found good agreement for several j, m! with the masses 
Eq. (2.18). We can also obtain the zero momentum bulk 



wavefunctions, from which we measure the dimensions of 
the dual operators, which are 2 + |. 

The real mess is found in the last set of scalars, dual to 
fluctuations of the gauge field in the internal directions. 
The difficulty is that there are three different vector har- 
monics on the internal space C3 that generally mix under 
the remaining symmetries of the problem. There are a 
few exceptions, including the j '■ — multiplct, for which 
the harmonic analysis is easy. These fields will feature 
prominently for the rest of the paper, as it is turned on 
when we have a background density and magnetic field. 
We therefore write their functional form for reference 

m' = : A = dtp + cos 9d(j), 
m! = ±2 : A = e ±l ^ [±id6 + sin 6d<t>) . (2.20) 

These fields are dual to the SO(3) x triplet scalar mass 

operator (QQ — QQ, QQ, QQ)- This operator has dimen- 
sion 1 in the supersymmetric theory, so some care needs 
to be taken in the bulk to ensure that we impose super- 
symmetric boundary conditions: the normalizable term 
of the near-boundary expansion behaves as a source and 
the non-normalizable is related to the expectation values 
of the dual operator. Unfortunately, we were not able to 
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find the spectrum of these, nor the other internal gauge 
fields. 

The meson masses we have obtained, while only rep- 
resenting half of the full spectrum, exhibit a huge de- 
generacy: all states of the same j + \m'\ + 4n have the 
same mass. This result is reminiscent of a broken 5*0(5) 
degeneracy in the D3/D7 system [S5], where all the me- 
son states were furnished by representations of 50(5) but 
not all states in an irreducible representation were degen- 
erate. We are finding hints of a similar effect here: the 
induced metric on the D6 branes is conformal to E 2,1 x C4, 

P[g] = (p 2 + m 2 )g 2 . x + 2 * 2 (dp 2 +4p 2 gc 3 ), (2.21) 

=9C 4 

where the C4 is a squashed, orbifolded 4-dimensional 
space. The conformal factor then depends on a coor- 
dinate on C4. Near the boundary the squashing vanishes 
and C4 also has local 5*0(5) isometry. We speculate that 
the hints of broken 50(5) invariance are arising from the 
global isometries of C4 near the boundary. The combi- 
nation of the conformal factor and the squashing pre- 
sumably breaks them to 50(3)r x U(1) x . It would be 
interesting to obtain the remaining meson spectra and 
more precisely understand what is happening on C4. 

3. Magnetic field and density 

We can now study chiral symmetry breaking in this 
theory. We will probe the phase diagram of this theory 
with a baryon density and magnetic field. The baryon 
current is dual to the U(l) gauge field living on the brane; 
we therefore add a baryon density and magnetic field to 
the field theory by turning on a field strength on the 
brane 

t->3 

F= 4k dp Adx ° + B dxl A dx2 ) ' ( 2 - 22 ) 

where we have normalized the field strength in the same 
units as the background metric Eq. |2.2| The field A is 



dual to the density operator and B is the magnetic field 
on the brane and in the field theory. We will consider 
brane embeddings that arc translationally invariant in 
the field theory directions (012) while only turning on the 
field strength above and the transverse scalar y = y(p). 
The equivalent ansatze in the field theory is that the 
correct extrema in the chirally broken phase are transla- 
tionally invariant and singlets under SO(3)r symmetry. 

These ansatzy however, are not consistent. There is a 
Chern-Simons term on the brane 

Scs = y / P[Ci] A F A F A F, (2.23) 

where C\ is the potential for 2- form flux F 2 . The pullback 
of C\ to our branes is 

P[d] = -k — . V (dip + cos OdiA (2.24) 

2vV + y 2 

so that this Chern-Simons term is linear in a flux on the 
internal space C3 . Careful study of this term shows that 
a density and magnetic field induces the (0; 1, 0) mode 
of the internal gauge field when the embedding is non- 
trivial. A consistent Ansatz therefore has a field strength 

4fcF 

— — = A / a (p)dpAdx +Bdx 1 Adx 2 +d(A tl) (r)(dip+cos6d(f))). 

(2.25) 

Recall that this internal gauge field is dual to the scalar 
mass operator. This field and y transform in the same 
representation of the global symmetry group and so the 
dual operators are both good order parameters for chiral 
symmetry breaking. It is interesting that both of them 
are turned on in this problem. 

Since our embeddings will have internal gauge fields 
turned on the brane action 

5 6 = -T 6 J d 7 He~^~P{g]+F + T 6 ]T / P[Q] A e F 

(2.26) 

will have three non-trivial Chern-Simons terms, namely 
those with C±, O3, and O7. For our ansatze the brane 
action becomes 



R9 

5 6 = -T 6 ^vol[IRP 3 ]vol[E 2 ' 1 ] 

v v ' 



dp 



pr x ll + y>i-A' 2 + -A' 2 X ll 



B 2 



1 + 4- 



y 2 (p + yy') _ vKba^ 

2r r 



3 A , A' 



(2.27) 



r 



where we have used 



P[0 3 



8 



r 3 dx° A dx 1 A dx 2 



(2.28) 



and have defined the positive orientation 
(a; , x , x 2 , p, ip, 9, 4>). For convenience, from here 
on we will consider the action density 56/volpE 2 ' 1 ] and 



give it the same name S§. Now we note that the field 
Aq only appears in the action Eq. (2.27) through radial 
derivatives. We therefore eliminate it in favor of a 



7 



constant of the motion, 



dC 

9l) 



= AfA' 



pr 



A 2 



1 + y' 2 - A' 2 



r± a/2 

p2 ^ 



yBA^ 



= d. 



(2.29) 

which we show in Sec. |II C"l] is the baryon density of the 
dual theory. In order to find solutions at fixed density, we 
algebraically solve Eq. ( 2.29 ) for A' Q and then Legendre- 
transform the action Eq. (2.27) to fixed density, 



S 6 = -Af dp 



- 4 ) + (dr - BA^y) 2 - 



y 2 (p + yy') 
2r 



3 A , A' 



(2.30) 



where d = d/N is a rescaled density. Solutions that ex- 
tremize this (fairly nasty) action will correspond to ex- 
trcma in the canonical ensemble of the dual theory. 



B. Little string theory 

The geometry supported a large stack of N D5-branes 
in type IIB string theory has a near-horizon limit |71] 

r R, 7 2 2 \ 

9 = + - {dr + r g §3 ), 

R r 



*<t> - 



9s R' R = ^ 9sN ' 



(2.31) 



where g s is the string coupling at infinity before taking 
the near-horizon limit and we have set d = 1. There 
are also N units of 3-form field strength on the 3-sphere 
sourced by the 5-branes. Type IIB string theory on 
this background is dual to the theory on the 5-branes, 
which in the g s — > limit is the (1,1) little string the- 
ory on a stack of N NS5-branes [71]. This theory is 
not a local quantum field theory and rather has a num- 
ber of stringy features, including string-like excitations, 
T-duality, and a Hagedorn spectrum at high tempera- 
tures [47] . For these reasons and because it lives in six 
rather than ten dimensions, it is called a "little string the- 
ory." For obvious reasons, the dual geometry is termed 
a "linear dilaton" background. The little string theory is 
strongly coupled in the ultraviolet and infrared free; in 
fact, the infrared theory is the maximally supersymmct- 
ric (5+l)-dimensional super- Yang-Mills theory on the D5 
branes. The full theory can also be obtained via decon- 
struction [72] . 

The linear dilaton background is conformal to E 6,1 xS 3 , 
i.e. to flat space. As a result, there is some difficulty 
in building a holographic map from bulk quantities to 
field theory observable^ The 5-brane background is 



rather special in this respect: the near-horizon region 
of the p-brane backgrounds are generically conformal to 
AdSp+2 x § 8_p , which yields a dictionary through holo- 
graphic renormalization |751 176j . 

Nevertheless we can learn a great deal about the lit- 
tle string theory from bulk physics. Bulk thermodynam- 
ics correspond to field theory thermodynamics; two-point 
functions in the field theory can be obtained by scatter- 
ing wave packets off of the center of the geometry [771 EE] • 
In this work, we revisit an item on this list: we can add 
flavor to the field theory by embedding probe branes in 
the dual geometry [5^1 E2] . 

There are a number of supersymmetric ways to add 
flavor to the little string theory, all of which are re- 
alized with dual probe branes that wrap cycles with 
four Neumann-Dirichlet directions with respect to the 
5-branes. These brane setups and thus the dual field 
theories preserve eight supercharges with flavor living on 
various codimcnsion defects in the field theory. The ex- 
ample we will look at in this work is that of codimension-2 
flavor. The weak-coupling, small N brane embedding has 
N D5 branes along the 012345 directions and flavor D5 
branes along the 012367 directions. In the large N limit 
the flavor branes wrap a (5 + 1) dimensional cycle inside 
the linear dilaton background including a circle inside the 
transverse 3-sphere. We write the lOd metric as 



ir 



R 



-95,i- 



- : {dp 2 +dy 2 +p 2 d9 2 +y 2 d0 2 ), 



R ^ffP + y* 

(2.32) 

where the radia l coordin ate r is related to these coordi- 
nates by r — yj p 2 + y 2 , the supersymmetric brane em- 



beddings are specified by x 



4 _ „5 _ 



<j) = constant and 



y = m, where m is proportional to the mass of the dual 
flavor. As in the ABJM theory, the 5-branes at the bot- 
tom of the geometry are located at p — y = 0. The 



2 There has been some work renormalizing asymptotically linear 



dilaton backgrounds from type II theories in 73 74 . 



boundary of the geometry is found as p or y —> oo. The 
field theory has a U(l) x U(l) R-symmetry, realized in 
the bulk as the U(l) x U(l) isometries of the two circles 
9 and <fi. Both symmetries are chiral. 

The field theory has an additional U(l) baryon num- 
ber symmetry which rotates the fundamental scalars and 
fermions. The baryon symmetry current is dual to a U(l) 
gauge field living on the flavor branes. Denoting its field 
strength as F, the action of the Nf flavor branes is given 
by 



S 5 = -NfT 5 / ^e-*y/-P\g\+F + S t 



cs. 



(2.33) 



where the branes have a tension T5 , the £ are coordinates 
on the worldvolume, P[g] is the induced metric on the 
branes, and Scs represents the couplings of these branes 
to the Ramond-Ramond form fields. 

As before, we will probe the phase diagram of this 
theory at nonzero baryon density and magnetic field. We 
therefore turn on a field strength 



F = A' (p)dp A dt + B dx 1 A dx 2 . 



(2.34) 



We will consider brane embeddings that preserve tran- 
sitional invariance in the field theory directions (0123), 
are flat in the defect directions (45), and simply wrap 
the ^-circle with <j> =constant. The only freedom left in 
these ansatze is to let the transverse position y of the 
branes depend on the radial coordinate p as y = y(p). 
The induced metric on the branes is 



P[g] 



vV + v 



R 



R 33,1 Vp 2 + y 2 

(2.35) 

The Chern-Simons terms with these ansatze are cubic in 
worldvolume fields that we do not turn on, and so we can 
consistently ignore them. Rescaling the spatial and radial 
coordinates by factors of R the brane action becomes 

S 5 = -NfT^vollS 1 } volpE 3 ' 1 ] 

S ✓ 



- : (dp 2 (l + y' 2 ) + p 2 d9 2 ). 



x j dp Pv /lT^T^y 1 + ^.(2.36) 



We pause to note that, at zero magnetic field, this is the 
same action for the D3/D3 system at nonzero density. 
Much of what we do (excepting the magnetic field) can 
then be recast in terms of that system. Also, as in our 
discussion in ABJM we will hereafter write of the action 
density S^/volpE 3,1 ] and give it the same name S5. 



Since the field Aq only appears in the action Eq. (2.36) 



through derivatives, we eliminate it in favor of a constant 
of the motion. Ordinarily this constant is interpreted as 
the baryon density [32l[80]. In this case, however, the 
gauge field has a near-boundary falloff A = plogr + 
0(1), where p is naturally interpreted as the chemical 
potential in the field theory. For this system (as with the 



D3/D3 system), the constant of the motion is instead 
proportional to p, 



dC 



B- 



dA' {p) ^1 + y'2 - A' 2 



Afp. 



(2.37) 



Solving Eq. (2.37) for A' Q and Legendre-transforming the 
bulk action to fixed p, we find the new action 



S 5 = -AfJ dpy/l + y' 2 




52 



(2.38) 



Curiously the bulk action at fixed p Eq. (2.38) corre- 



sponds to the field theory action at fixed density. We 
will return to this in Sec. Ill C 21 

There is one more control parameter that we can triv- 
ially turn on: the analogue of the scalar mass we looked 
at in the ABJM theory Eq. ( |2.20| . In the D3/D3 system 
as well as here, the dual field is the zero-mode of the in- 
ternal gauge field on the ^-circle. Turning it on with a 
field strength 

F = A' Q (p)dp A dt + B dx 1 A dx 2 + A' e {p)dp A d6, (2.39) 
the brane action becomes 



S 5 = -M J dppJl + yii 



A' 2 



~A'h 1 



B 2 



(2.40) 

As before, we note that Ao and Ag both appear through 
derivatives so that we have two constants of the motion. 
The reader can verify that Ag has the large-p falloff Ag = 
M log p + 0(1), so that the constants are 



dC 



AfpA> J 



1 



9A' (p) 
dC 



^l + y^-A' 2 + ^A' e 



= Afp, 



Nr 2 A> B J 



1 + 



B- 



pJl + y'2-A'i 



1 A' 2 



: MM{2A\) 



Solving for A' and A' e in terms of p and M we find 
p^l + y 12 



A' 

4' 

R 2 



f'- 



il 



B 2 -M 2 - 



Mp 2 y /l + y r - 



B 2 -M 2 " 



(2.42) 



Legendre transforming the bulk action Eq. ( |2.40 ) with 
respect to both A' Q and A' e we find the new action 



S 5 = -Mj dpVl + y' 2 Jp 2 + p 2 (l 



B 2 -M 2 ' 



(2.43) 
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Two comments are in order: (i.) the effect of adding 
a scalar mass can be completely solved in terms of a 
magnetic field with mag nitudc B = VB 2 - M 2 , 



new 



(ii.) reality of the action informs us that there is an upper 
bound on the scalar mass, M^ ax = p 2 + B 2 , that can be 
described with probe branes. 



C. Holographic renormalization 

1. Probe branes in ABJM 

We continue by holographically renormalizing flavored 
ABJM theory. We have several tasks in this subsection: 
(i.) to show that the embeddings we propose are dual to 
massive flavor have zero free energy, (ii.) to obtain one- 
point functions at nonzero density and magnetic field, 
and (hi.) determine the proper operator normalization 
for the transverse scalars. We will achieve all of these 
tasks via holographic renormalization. 

The action Eq. (2.30) diverges near the AdS4 bound- 



ary. This infrared divergence is related to an ultraviolet 
divergence of the dual theory. We consistently renormal- 
ize it by diffeomorphism-invariantly regulating the bulk 
theory. We do this by first solving the equations of mo- 
tion for y and near the boundary, 



n=l y 



A n 



(2.44) 



where the higher y n are recursively determined by m and 
yi and the higher A n are determined by Ai and A^. Next 
we integrate the action on an arbitrary solution up to a 
very large cutoff p = A, 



<S 6 ,a = -JV / dpp 2 + finite. 



(2.45) 



The brane action then has a single divergence from the 
volume of AdS4. We regulate it by adding a local coun- 
terterm on the cutoff slice 



5, 



CT 



lp=A, 



(2.46) 



where 7 is the induced metric on the cutoff slice. For an 
action with a more general divergence pattern we would 
add a number of local, (cutoff-)diffeomorphism-invariant 
counterterms on the cutoff slice. 

We define our renormalized action by adding the coun- 
terterm to the action Eq. (2.30) and then taking the cut- 
off to infinity, 



lim 

A— ¥00 



S: 



5,A 



'CT 



(2.47) 



This bulk action is finite and convergent on-shell; from it 
we can obtain correlators of the dual theory in the usual 
way. We begin with our first goal, the free energy of the 



embeddings y 
branes is just 



to. The renormalized action for these 



'6, rcn 



= -AT- 



(i 



(2.48) 



which is nonzero. However, we have missed one subtlety: 
we have not specified the correct infrared boundary con- 
dition for C7. To see this, recall that there is a term 
in the brane action J P[C 7 ] that changes under gauge 
transformations C 7 n- C 7 + dC 6 with the boundary term 
S n- S + J d Cq , where d is the boundary of the worldvol- 
ume. Gauge transformations with support at the world- 
volume boundary change the boundary conditions on the 
bulk theory. In particular, when Cq does not vanish on 
the AdS4 boundary it changes the source for a magnetic 
current in the dual theory [ST]. On the other hand, a 
non-trivial Cq at the bottom of the brane changes the 
infrared value of the gauge field. We will impose the 
boundary condition that the gauge field vanishes there. 
However, the pullback of C 7 in Eq. (2.15 1 does not vanish 
at the bottom of the brane, 

R 9 

P[C 7 ]| P=0 = y^m 2 sm6dx° Adx 1 Adx 2 AdrAdipAdeAdcf). 

(2.49) 

In order to be consistent with our boundary conditions 
we therefore need to perform a gauge transformation so 
that P[CV] vanishes at the bottom of the brane without 
changing its near-boundary behavior. This can be done 
with a Cg that vanishes at large r and at smaller r be- 
comes 



Ca 



7? 9 



2 8 fc 2 3 



sin 9 dx° A dx 1 A dx 2 A di/; A d6 A d(j). (2.50) 



This transformation adds a contribution to the bulk ac- 
tion proportional to m 3 . In fact, the renormalized action 
in this gauge simply vanishes. The free energy —S of the 
dual background thus vanishes as advertized. 

Next we will compute one-point functions. In order 
to obtain the density, let us consider the brane action 
before Legendre transforming Eq. (2.27). This action 



diverges in the same way as its Legendre transform and 
so we define S^rcn i n the same way as before. This action 
is equal to (minus) the free energy of the dual state in 
the grand canonical ensemble. Under a small variation 
A Q h-> Ao+SAq, y h> y+Sy, A$ H> A^+SA^, the on-shell 
action changes by a boundary term on the cutoff slice 



_ _ . dC . dC dC 
S hdy = SA 0M r + 6y— + SA il -^r. 

Recalling the near-boundary expansion of A 

d 

A = »- Jrp + ..., 



(2.51) 



(2.52) 



and treating the leading term A\ in the large-p solution 
of A^p as the source for the dual operator, we therefore 
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find the one-point functions 



= (Oy) = 



bdy 



dm 



d = <i°) = 
C = (0 4 ,) = 



dSbdy 

dp 

dSbdy 

~ dAi 



= -AfA 2 



(2.53) 
(2.54) 
(2.55) 



However as we noted in Sec. |II A2[ these boundary con- 
ditions on A^p do not correspond to the supersymmetric 
boundary conditions where A 2 is the source for the dual 
operator. To implement this boundary condition we need 
to add a total derivative to the bulk action 



^6,1'cn 



i — ^ s, 



6,rcn 



^ J dp ( P * A%)' 



(2.56) 



which accordingly shifts the boundary action. The one- 
point function of O^p with these boundary conditions is 
then 



C={0 4 ) = - d ^=NA 1 . 



(2.57) 



Notably, all one-point functions vanish for the embed- 
dings y — m as advertized earlier. 

Finally, we will obtain the relative normalization be- 
tween the hypermultiplet mass operator M. and O y . We 
could follow |82j and compute the two-point function 
(Oy(x)O y (0)) from holographic renormalization as well 
as (M(x)M(0)) from the field theory. We elect to take 
a simpler approach and compute the mass of a straight 
string that hangs from the bottom of the flavor brane to 
the D2 branes at the bottom of AdS2, 



1 f m 

M = — / dp^P[g~] = 
27T Jo 




(2.58) 



This is the mass of the lightest charged quasiparticle in 
the dual theory and therefore corresponds exactly to the 

hypermultiplet mass. We therefore have A4 = \J^O y . 
This relative normalization contains an important piece 
of physics: the quark mass (M) is larger than than meson 
masses (set by the scale m) by a factor of ^N/k, the 
square root of the 't Hooft coupling of the ABJM theory. 
This is precisely analogous to the mesons in the D3/D7 
system (among others), where the meson masses go like 
«^quark/"\/A, where A is the 't Hooft coupling of the gauge 
theory [52] . 



2. Probe branes in the linear dilaton background 

Two comments are in order before we continue. First, 
because the linear dilaton background is conformal to flat 
space rather than to AdS, the regularization scheme we 
propose should be considered to be just that: a regu- 
larization scheme. We have not rigorously renormalized 



the bulk theory. As a result, our main reason for study- 
ing this regularization is to nail down the proper holo- 
graphic renormalization in the D3/D3 system at nonzero 
density [55]. The technical details of that problem are 
identical to those we face in creating our regularization 
scheme. 

While our regularization and so our free energy for 
the dual theory has issues, our prescription for the chiral 
symmetry breaking condensate is more reasonable. Our 
result for this condensate is identical to its computation 
in the D3/D3 system and simply returns the normalizable 
term in the large-p expansion of the field y. This matches 
the usual AdS/CFT picture where the normalizable term 
of a field is proportional to the vacuum expectation value 
of the operator sourced by its non-normalizable term. 



An embedding that minimizes the action Eq. (2.43) 
will have the large-p falloff y = m log p + yo + 
0(log 2 p/p 2 )- Integrating the action Eq. (2.431 up to a 
large radial cutoff A 



find 



S: 



5.A 



Af dpp(l + 



B 2 -M 2 



2p 2 



■ finite. 



(2.59) 

The brane action then has a volume divergence from the 
first term and four logarithmic divergences from the four 
worldvolume fields we have turned on. We can regulate 
these divergences with local counterterms on the cutoff 
slice p — A in a similar manner to holographic renormal- 
ization, 



(2.60) 



with the counterterms 

A = 1 
c 2 = 

Co = 



£4 = 



u = - 




An 



(2.61) 



4 

1 

2 log/ 

where 7 is the induced metric on the cutoff slice and all 
induces are slice indices contracted with 7. Note that in 
the D3/D3 system the proper counterterm for the gauge 
field £2 does not have a factor of the dilaton. We then 
define the regulated action 



S^rcn = lim 

A— >oo 



£5, A + SCT 



(2.62) 



We should contrast this result with the divergences of 
the original action Eq. (2.401, 



>5,A 



= -Af dppll 



m 2 ~ p 2 + B 2 
2p 2 



M 



(2.63) 
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which differs from the diverg ence s tructure of the Legen- 
dre transformed action Eq. (2.43) in the sign of the p 2 
and M 2 logarithmic divergences. The Legendre trans- 
form changes the action's asymptotics and so we need 
slightly different slice countcrtcrms, 



with 



C 4 



£i = £i, £2 — £2, £4 
£3 = —£3, £5 = —£5, 
from which we define the regularized action 

-S^ren = lim i S 5,A + Sct] ■ 
A— foo 



(2.64) 



(2.65) 



(2.66) 



What physics is contained in these divergences? First, 
the logarithmic divergences of the bulk action correspond 
to a Weyl anomaly in the dual theory proportional to 
terms like m 2 , p 2 , &c. We can see this by performing 
an infinitesimal Weyl variation in the boundary theory, 
which corresponds to an infinitesimal scale transforma- 
tion in the bulk. The logarithmic terms in the bulk action 
transform inhomogeneously, giving rise to a change in the 
action under the variation [S3]. Second, some terms in 
the anomaly differ by a sign depending on the field theory 
ensemble. 

In order to make this statement more precise we need 
to relate bulk fields to expectation values in order to see 
what field theory ensembles correspond to our regular- 
ized bulk actions. Usually the original action Eq. (2.40) 



corresponds to the grand canonical ensemble and the 
Legendre-transformed action to the canonical. In order 
to see if this is the case, we need a prescription for com- 
puting expectation values. 

Even though we do not performed strict holographic 
renormalization for the bulk theory, let us proceed and 
define one point functions in the field theory by variation 
of our regularized actions. Moreover we will assume that 
the original action Eq. (2.40) corresponds to the grand 



canonical ensemble and then show that this assumption 
is thermodynamically consistent. Varying the original 



action Eq. (2.66) with respect to the fields y, Aq, and Ag 



we find a variation of the boundary action on the cutoff 
slice 

SS hdy = Af[ P {-Sy y' + SA A' - 6A e A' g ) 

+^ x (5yy-5A A + SA e A e )}. (2.67) 

We define our one-point functions by variation of this 
boundary action with respect to the sources m, p, and 
M in the A — ► 00 limit, for example 



(Oy 



SS ] 



bdy 



,. , . <5<Sbdy 

lim log A — - — -. 

A->oo oy 



(2.68) 



Recalling the large-p expansion of our worldvolume fields, 

y = mlogp + yo + ^0 = plogp + a + 
A 9 = Mlogp + co + ... (2.69) 



we therefore find our one-point functions 

SSbdy 



C = (Oy) 

C = (Og) 



Sm 

SSbdy 

Sp 

SShdy 

SM ' 



: A/C . 



(2.70) 



What happens when we Legendre transform? The 
transformed action Eq. (2.43) is related to the original 
Eq. ([2~40| by 



5*5 — 05 



Ag 



d£ 
OA' 



(2.71) 



Since each of the transforms gives a logarithmic diver- 
gence at large p the counterterms required to regularize 
S5 are different than for S$. Noting that ag and Cq can 
be obtained by integrating the gauge fields, 



a 



lim 

A— >oo 



c = lim 

A— >oo 



' .A 

/ dpA' Q -p\ogA 
■ ,A 

/ dp A' g - M log A 



(2.72) 



the reader can show that S$ jren is related to S^ren by 



S. 



5,ren 



SI 



5,rcn 



pd - MC. 



(2.73) 



Since the bulk action corresponds to (minus) the free 
energy of the dual equilibrium state, we have 



F = F + pd + MC, 



(2.74) 



so that the original bulk action S$ indeed corresponds to 
the grand canonical ensemble and S$ to the canonical. 

One nice check for these results is to compute the den- 
sity in the grand canonical ensemble for the symmetric 
phase c = 0. We can do this in two different ways: (i.) 
by our prescription for the density in Eq. (2.70) and (ii.) 
by computing the free energy F = —Se. icn for this phase. 
Thermodynamic consistency of our regularized free en- 
ergy requires that d = —dF/dp. The regularized action 
and dual density for the symmetric embedding y = can 
be evaluated by means of Eqs. (2.66 )(2.70), and (2.72) 
to be 



^6,ren 



r h Jp 2 +B 2 



-(p 2 + B 2 



1 - 2 log 



p 2 + B 2 + r 2 



d =plog 



B 2 



(2.75) 



The reader can verify that d — dSe tren /dp, so that our 
prescription for the free energy is indeed thermodynam- 
ically consistent in the symmetric phase, at least with 
respect to density. 
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III. RESULTS 



ABJM with flavor 



Our goal in this subsection is to identify the existence 
and location of the BKT transition in this system. As 
a result, we will not solve the full problem of obtaining 
embeddings, condensates, &c in the broken phase. The 
full problem, while soluble, is numerically daunting for 
little payoff: the final results will look very similar to 
those of the D3/D5 system, but the broken phase will be 
described by two condensates rather than one. 

As in the D3/D5 system, the onset of the chiral BKT 
transition can be found by analyzing the stability of fluc- 
tuations around the symmetric embedding y = = 0. 
The se fluc tuations are governed by the quadratic part of 
Eq. (2.301, which is equivalent tq^l 



£ 6 ~ -y [Y' ■ K ■ Y' + Y ■ V ■ Y] 



where 



Y 



V 



( y/p 

I A 



,2 P 4 +B 2 
v /= G 
dB 
v /3 G 



-G P 2 I 2 , 



dB 



-G 



-3p 



2 , £ 



- 4 + B 2 ' 



f—G 



and 



d 2 + B 2 



(3.1) 

(3.2) 
(3.3) 

(3.4) 



is the determinant of the open string metric on the sym- 
metric embedding. This system behaves rather differ- 
ently in two different regimes. At large p ^> y/B, the 
fields y/p and decouple from each other and fluc- 
tuate as m? = —2 scalars in AdS2- Fluctuations here 
correspond to the ultraviolet of the dual theory, while 
fluctuations at small p <C \J~B~ correspond to the deep 
infrared. Here y/p and A^, fluctuate as coupled scalars 
in AdS2- The infrared eigenmodes are 



$± = A 4 , + 



B±\fd 2 



JPy 
P 



(3.5) 



with masses m± = ^B/ y d 2 + B 2 in the effective AdS2 
region. This is the crux of the argument: while the D6 
brane is always embedded in an AdS4 background, the 
fields living on the brane see different geometries at dif- 
ferent p. The open string metric [84] on the brane is 
dominated in the infrared by the control parameters and 
there becomes AdS2-like (the full result is a Lifshitz-likc 
space, but this does not affect our results). 



3 By equivalent, we mean that there are terms of the form f(p)yy' 
in the quadratic action whose contribution to the equation of 
motion can be compensated for by a term — - y 2 . 



For large density d 2 / B 2 > 15, the mode $ S gn(B) fluc- 
tuates stably in the IR. However for d 2 / B 2 < 15, the 
mass of 4 > sgn (B) passes below the Breitenlohner-Freedman 
bound for stability ^49 in the effective AdS2, TObf = 
— 1/4. There will be an infinite, exponential hierarchy 
of tachyons in the symmetric embedding signalling an 
instability toward a chiral symmetry broken y, A^, =/= 
embedding. There is therefore a chiral quantum phase 
transition at 



d 

W\ 



N 



2 5 / 2 7T 



15. 



(3.6) 



The appearance of the factor N N/k is analogous to the 
factor of Ny/X that appears in the chiral transition in the 
D3/D7 system [21], where A is the 't Hooft coupling of 
that gauge theory. 

We have verified the existence and location of the chi- 
ral quantum transition. On the grounds of [22, 50 , we 
expect that it is of the holographic BKT type. We can 
back up this claim with more rigor by sketching a deriva- 
tion of the stable embeddings in the broken phase. These 
embeddings determine the thermodynamics of the bro- 
ken phase. Near the transition we find the embeddings 
by matching solutions in two regions: a small-p non- 
linear core and a large-p linear bulk. We begin with 
the core region p <C V~B where the brane action is well- 
approximated by 



5 6 = -M jdp-Jl 



,/2 



— A' 2 



p 2 B 2 {\ 



(dr- BA^y) 2 .(3.7) 



Our brane embeddings extremize this action and obey 
the infrared boundary conditions y(0) = with A^, 
falling off so that the action is integrable near p — 0. 
Importantly, this region enjoys a scaling symmetry 

(3.8) 



As a result we can find a family of embeddings that obey 
our boundary condition: for y — f y (p),A^ = /^(p) a 
particular solution that obeys our boundary conditions, 
we define the solutions y^(p) = £,f y {p/Q,A^£ = U(p/0 
for all £ € R + . However, this scaling symmetry is not 
enough to obtain the two-parameter family of core solu- 
tions. We need an extra boundary condition, which can 
only be imposed by matching to the near-AdS4-boundary 
asymptotics. Moving on, consider large p in this region 
(near the boundary of the IR AdSj), y and A^ mix in the 
same way as we noted in Eq. ( |3.5[ ); there are two scalars 

with the asymptotics of m 2 = ±B /yd? + B 2 scalars in 
AdS2- We are interested in embeddings close to the chiral 
transition, so we consider \B\ close to the critical value. 
Defining 



a± = 2 



'BF 



= J-l± 



AB 



VB 2 + d 2 



(3.9) 



13 



this region corresponds to a asn m) both real and small. 
Now letting B > 0, for a particular solution y,A^, the 
eigenmodes Q± have the near-AdS2-boundary asymp- 
totics 



$+0) ~ p - 



-a cos (| log pj 



2ai 



(flogp 



(3.10) 



where the a.; and 6^ depend on a = a + and we have 
defined a = ia. There will be a one-parameter family 
of ai's and frj's that, together with scaling, generate the 
asymptotics of all core solutions. 

To order a, a general embedding that obeys our bound- 
ary conditions will have the large-p asymptotics 



^+(p) = ^sin(flog^ 
V p a \ 2 pi 



, p > pi = £e 



_ p aa/ax 



$cAp) = y ^ ( & o£ f P 5 + &i£ 5 p 5 ), P»Pi (3.H) 

The can be series expanded in a 2 as 



«o(a) = X! a o," a2 ": 



n=0 



f'l 



ai, n a 



2/i 



(3.12) 



n=0 



In principle, we could numerically measure the aj iTl by 
obtaining core solutions that obey our boundary condi- 
tions. However, we elect to continue and leave our results 
in terms of the aj, n . 

The linear region is the domain where y/p,y' , A^, and 



pA', -C 1 and so the quadratic Lagrangian Eq. (3.1) is 



a good approximation to the full Lagrangian. Unlortu- 
nately, we have not been able to decouple the fields y and 
Aty everywhere in this regime and thus have not found the 
linearly independent solutions for them. However, we do 
know the solutions for them for p <C y/~B in the effective 
AdS2- The eigenmodes $± fluctuate there as scalars of 

mass-squared ^fB/y dP + B 2 . The effective AdS2 region 
extends from p of order but smaller than s/B to where 
our linearized approximations begin to fail. 

For vanishing sources and small condensates we can 
further approximate the solutions in the linear region: 
we linearize in c and C. Near the transition the solutions 
in the effective AdS2 region will then be 



$ + (p) ~ p 
$_(p) - p-3 (/ 3 pf + / 4 p"f ) , 



- sin (- log pj + / 2 cos (- logp ) ) , 



(3.13) 



where the fi are linear in c and C and depend upon a. 

The linear and core regimes overlap in the region p <C 
VS, $±,y', and pA'^ <C 1. This region is rather large 
near the transition and includes the infrared AdS2 region. 
We match the linear and core solutions at an arbitrary 



point in this region, which we take to be of order £, is 
near the bottom of the AdS2- Both the core the linear 
asymptotics ov erlap here so that we need only match the 
solutions Eqs. (3.11) and (3.13). In order for these to 



match the argument of the sin in Eq. (3.13) must make 
7r between p\ ~ £ and p ~ y/B. To exponential accuracy 
we then have £ ~ e~ 27r / a which gives 



(3.14) 



Thus the chiral transition has exponential scaling and is 
of the holographic BKT type. 

Next, the way that we matched the core and linear 
solutions makes it clear there are in fact an infinite num- 
ber of symmetry-breaking embeddings. Following [22] we 
term these solutions "Efimov extrema" for their resem- 
blance to Efimov states [S5]. The most general match 
will have the argument of the sin in Eq. (3.13) can pass 
through row between p\ ~ £ and p ~ vB for n any posi- 
tive integer. These embeddings oscillate through n/2 pe- 
riods between the D2 branes at the bottom of AdS4 and 
the boundary. Matching the core and linear solutions at 



p ~ £ we find to exponential accuracy 
which yields condensates 



t/c 



-2mr / a 



(3.15) 



where C\ , C\ are the condensates for the simplest em- 
bedding above. These embeddings, however, are not 
dual to ground states: denoting the free energy of the 
symmetric phase as Fq, we find that the free energy for 
the state dual to the n th non-trivial embedding goes like 
F -F n - e - 27m/a . It then follows that the n = 1 embed- 
ding is dual to the ground state in the broken phase. The 
extra non-trivial embeddings are then dual to extrema of 
the dual theory that are at best metastable. The free en- 
ergies of these extrema form an infinite tower, spaced by 
the same factor e -27r / Q , reminiscent of the Efimov states. 
These are an infinite exponential hierarchy of three-body 
bound states in atomic systems where two-body interac- 
tions have been tuned to threshold. 

Finally, some words about finite temperature are in 
order. We turn on a temperature in the field theory by 
heating up the D2 branes. The dual geometry is a black 
brane geometry with an AdS4-Schwarschild factor and a 
horizon radius r/, oc T . This change will dramatically 
affect the chiral transition. In the D3/D5 system the 
chiral transition became second-order with mean-field ex- 
ponents at any nonzero temperature [22, 45 . Indeed, we 
expect the same thing to happen here. 

The crux is that fluctuations on the brane no longer 
see an infrared AdS2 geometry arbitrarily close to the 
horizon. Rather, the equation of motion for fluctuations 
will have three separate distinct regimes: an AdS4 region 
at large rlVB, an infrared AdS2 for Th -C r VB, and 
a near-horizon geometry conformal to Rindler space for 
r — r/i <C Th- The fields $± are therefore stabilized in the 
deep infrared. Their fluctuations can be described by 
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an equivalent quantum mechanics problem with a non- 
trivial potential [50]. This potential has the form of an 
inverse square potential regulated at both short distances 
(by the control parameters) and at large distances (by the 
temperature). As d is decreased, the coefficient of the 
inverse square potential passes below the stable bound; 
however, no tachyons are formed because the potential is 
regulated in the IR and UV. However, as d is decreased 
further a single tachyon will eventually form signalling 
a second-order continuous transition rather than BKT 
scaling. This is very different from the picture at zero 
temperature: there is no infrared regulator there and so 
an infinite exponential hierarchy of tachyons forms at the 
transition. Thus the holographic BKT transition only 
occurs at zero temperature, as generally noted in [221123] . 

The exponential scaling, however, is not utterly lost at 
nonzero temperature. As shown in Fig. [T] for the little 
string theory, the finite temperature condensates asymp- 
tote to the zero temperature result for large enough mag- 
netic fields. Moreover, the critical temperature will itself 
scale exponentially with 1 j a 23J . The simple way to see 
this is to use the quantum mechanics picture for fluctua- 



Eq. (2.431 



tions. For a scalar field in AdS^+i with m 2 



-a 2 /4 



with ultraviolet and infrared regulators that become im- 
portant at respective scales ruv and Hr, the first tachyon 
forms when they are exponentially separated [50 



a ruv 
2 r m 



(3.16) 



The AdS2 region is defomred in the infrared by the black 
brane horizon at ~ r^. To exponential accuracy the 
finite temperature transition therefore occurs at 



-2ir/a 



(3.17) 



These features should be quite general for holographic 
BKT transitions. 



B. Flavored little string theory 

This subsection is divided into three parts. In the first, 
we do some analytic work and identify the BKT transi- 
tion in this system. Next we obtain some numerical re- 
sults at zero temperature and find excellent agreement 
with our analytic predictions. Finally, we look at this 
theory at finite temperature wherein we lose the BKT 
scaling near the transition. Our discussion throughout 
will parallel our work above in the ABJM theory. 



N^fp 2 



B 2 



-II ' 



MB 2 y 2 



B 2 



P 

(3.18) 

We pause to consider the behavior of y in two different 
limits. At large p ^> p,B,M the radial equation of mo- 
tion for y/p is that of a m 2 = — 1 scalar in AdS3. At 
small p (the infrared of the dual theory), however, y/p 

behaves like a m 2 = — n B ~ scalar in AdS 2 . This re- 

suit is crucial: even though the D5 branes are embedded 
in the linear dilaton background at all p the fields liv- 
ing on them see different geometries at different p. The 
open string metric is dominated in the IR by the control 
parameters |84) . 

For p 2 /B 2 > 3, y/r fluctuates stably in the IR. How- 
ever, for p 2 /B 2 < 3 the mass of y/p in the IR is below 
the Breitenlohner-Freedman bound for stability [49] in 
the effective AdS 2 , m| F = —1/4. The y — embed- 
ding is unstable in this phase and the stable phase will 
be y 5^ 0. There is therefore a chiral quantum phase 
transition at 



EL 

B 2 < 



= 3. 



(3.19) 



We can compute thermodynamics in the broken phase 
by solving for the stable embeddings with appropriate 
boundary conditions. Near the transition, we can find 
the embeddings by matching solutions in two regions: a 
small-p non-linear core and a large-p linear bulk. We 
begin with the core region p <C p, B, M where the action 
looks like 



-N / dp^/l + 



y' 2 \ p 2 



B< 



(3.20) 



Our brane embeddings minimize this action and obey 
the boundary condition y(0) = 0. Moreover, since this 
region enjoys a scaling symmetry we can find all such 
embeddings: for f(p) a particular solution that obeys the 
boundary condition, we define another solution y^(p) = 
£/(p/0 for alU e H+. 

At large p (near the boundary of the IR AdS 2 ) y/r 
has the asymptotics of a m 2 = — J 3 ^^ scalar in AdS 2 . 
Since we are interested in the near-critical behavior, we 
will consider p 2 close to the critical value. Defining 



"BF IR 



p 2 + B 2 ' 



(3.21) 



1. Analytic results 



this region corresponds to a both real and small. Then 
/ has the large-p asymptotics 



The onset of the chiral BKT transition can be found 
by analyzing the stability of fluctuations of the symmet- 
ric embedding y — 0, given by the quadratic part of 



/{/>) -- v p { -a cos (I logp) + ^ sin (~logp 



(3.22) 
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where ao and a\ are functions of a. To order a, a general 
embedding then has the large-p asymptotics 



V(P) = —Pepsin ( £ log f-Y p > Pl = ee ao/Ql . 
a \2 Pi J 

(3.23) 

Since we are interested in the near-critical embeddings, 
we series expand Op and a\ in a, 



a (a) = J3 ao^a 2 ", 



n=0 



ai(a) = ^ «i. 

71 = 



2/i 



(3.24) 



As we discuss in Sec. |IHB"2| we found solutions / at var- 
ious small a and numerically measured ao.o = —.32189, 
a , 2 = -2.328, oi, = .42525, and oi, 2 = -.7477. 

The linear region is the domain where y/p, t/ < 1 and 

is a good approxi- 



the quadratic Lagrangian Eq. (3.18 



inearly independent 



mation to the full action. The two 
solutions for y can be found here without any further 
approximation 



i±i« f l±ia l±ia ia , 
/±(«)=« » 2 Fi ( — ^— ■l±—--u 2 



where 



p? + B 2 



(3.25) 



(3.26) 



Since we are studying the theory at zero bare mass, we 
impose the large— p boundary condition 



y(u) = c+f+(u) + c_/_(u) =y + O 



(3.27) 



where we recall from Eq. (2.701 that the condensate of 
the dual theory is simply —Nyo- Demanding the reality 
of the embedding forces c± to be complex conjugates, 
which together with the boundary condition Eq. (3.27) 



determines 



2n 



c± = I Ti C 2n+iu 2n 1 + J3 C 2n+2 cr 

\ n=0 n=0 / 

(3.28) 

where the Ci are all real constants. The hrst four are 



Ci 



V2 C 2 



d 



log 8 Cg = 
C x 

1 fCi 
Ci 



2 

C2C3 



9i 
c x 

C(3) 



7T" 

24 ' 



(3.29) 



C? 3 VCi/ 4 
At small u near the IR AdS 2 boundary y takes the form 



y(u) =Vu{c+u ia / 2 + C-U- ia / 2 ), 



(3.30) 



which can be rewritten with Eq. (|3.28|) as 

2CiCy/u 1 



Afc 



E^2n+l 2n 
— — — or sin 

n=0 



n=0 



C 2 n- 



logu (3.31) 



2 „2n+i 



(flogu) 



To order a, these asymptotics can be rewritten as 
2C lC ^p 



y(p) 



Na(p 2 + B 2 ) 1 / 4 



where 



P2 



p 2 + B 2 exp 



a p 

2 P2 



2C 2 
' Ci 



(3.32) 



(3.33) 



is an intrinsic length scale in the linear regime. 

The linear and core regimes overlap in the region 
p <C Pi y/p,]/ <C 1. This region is quite large near 
the transition. We therefore perform the matching at 
an arbitrary point in this region, which we pick to be 
p = £ where the core and linear asymptotics overlap. In 
order to match these solutions the argument of the sin 
in Eq. (3.32) must make 7r between pi ~ £ and p 2 ~ p. 
To exponential accuracy we then have £ ~ e -2n/a wn \ ca 
gives a condensate 



-it J a 



(3.34) 



Thus the transition in the little string theory obeys BKT 
scaling. The exponent and normalization of both £ and 
c can be found by matching y and y' at £ order by order 
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£ = y p 2 + B 2 exp 



2tt ag 2C2 o 

+a 

a ai Gi 



V 



12af~ + 3Cf ^ 2 + ^i^ 2 ^ 3 _ 3CiC4) 



V 



+ 0(a 4 ) 



/ 



/ 



p 2 + B 2 ^- exp 



'a 2ai Ci 1 2 



no 



2C\ 



(Ci + 2dC 3 ) +0(a 4 ) 



(3.35) 



We also compare this result to numerical data in Fig. [T] 
As before, the way that we matched the core and linear 
solutions makes it clear there are an infinite number of 
embeddings which, following we again term "Efimov 
extrema" for their resemblance to Efimov states [55] . In- 
deed, the most general match has the argument of the sin 
in Eq. ( 3.32 ) pass through rnr between pi ~ £ and P2 ~ /U, 
for n a positive integer. These embeddings go through 
n/2 oscillations between the bottom of the brane and the 
boundary. Matching y and y' at p = £ we find a scale 
factor and condensate corresponding to these extrema of 



\n+i- ( "~, 1),r „ 1 /^/„4^ 



(-l) n +V 



■ci + 0(a 4 ), 



(3.36) 



where £1 and c\ are the scale factor and condensate for 
the simplest embedding above. However, the additional 
extrema are not ground states: if one denotes by Fq the 
free energy corresponding to the symmetric embedding, 
then the free energy of the n th non-trivial extremum goes 
like F — F n ~ e ~2mi/a^ SQ tne n — 1 embedding is 
the ground state. The infinite tower of extrema, spaced 
by the same factor e~ 27r / Q is reminiscent of the Efimov 
states. 



2. Numerical results 
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FIG. 2. The free energy difference between the broken 
and symmetric phases of the little string theory as a func- 
tion of magnetic field. This curve scales exponentially as 



AF ■ 



at the bottom of the brane. This corresponds to <f> = 4>{t) 
going to zero as r — > oo. Indeed, there is a one-parameter 
family of series solutions that obey this boundary condi- 
tion, 



-) = fae- mT 1 + E ' 



-2mr 



(3.39) 



We have two tasks in this subsection. The first is to 
identify how we numerically obtain core solutions (i.e. 
solutions that minimize Eq. (3.20)). The second is to 



numerically obtain full zero-temperature numerical so- 
lutions that minimize the whole action Eq. (2.43). We 



will obtain numerical solutions at finite temperature in 
Sec. UJJJJ3J 

We find numerical core solutions by shooting from 
small p, i.e. at the very bottom of the brane. In or- 
der to find a good series solution here we elect to change 
coordinates and redefine our worldvolume fields by 



V 



P : 



(3.37) 



The core action in these coordinates is 



'J dre' T ^f\ 



'^/p 2 + B 2 sin 2 (f>. (3.38) 



S 5 - -. V / ,W< r v 1 

As the reader may note from the form of the behavior of 



y in the AdS2 region Eq. (3.30), both y/p and y' blow up 



where 



-1 



|B 2 



(3.40) 



and the fa can be recursively solved for in terms of fa. 
Since this series expansion is in powers of e~ mT we could 
write a series solution for <f> in terms of the radial coordi- 
nate r as <f) = far m + 0(r 3m ). In any event, we use this 
small— p series solution as initial data upon which we nu- 
merically integrate to find y up to very large p close to the 
AdS2 boundary (of order p ~ 10 25 ). We then match the 
numerical solution to the near-AdS2 boundary solution 
Eq. (3.23) to find ag and a\ at small a. 



Our full zero-temperature solutions were obtained by 
shooting from large p. Imposing the zero mass condition, 
we find a series solution for y near infinity 



y(p) = Vo 



VnP 



-2n 



(3.41) 
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FIG. 3. The density as a function of chemical potential in 
the broken phase of the little string theory. The thin dashed 
blue curve represents a direct measurement of the density by 
Eq. (2.701 while the thick dashed black curve is —dF/dp. 



where the y n are recursively determined by yo. We use 
this series solution as initial data upon which we again 
numerically integrate. This time we proceed from large 
p (here, we use a series solution to eight orders and inte- 
grate from p = 100) to very small p (of order p ~ 10 
to 10~ 28 ). However, due to the spiky behavior of the core 
solution Eq. (3.39) near p = we need to more careful in 



.">o 



order to ensure that we correctly shoot to the IR bound- 
ary condition y(0) = 0. This simply means that we shoot 



for our solutions match the core series solution Eq. (3.39) 
at small p. 

We plotted some of these results in Sec. [T] First, wc 
plotted the condensate in the broken phase in Fig.[l] Wc 
also computed the regularized free energy F = —Ss >ien 
in the broken phase and plotted it in Fig. [2] The agree- 
ment between our analytical and numerical data for the 
condensate indeed holds to 0(a 4 ). 

Finally, we recall our prescription for computing the 
density d in the dual theory. In principle wc can also 
obtain the density from the free energy: thermodynamic 
consistency requires that d = —dF /dfi. Wc compare 
these in Fig. [3] after using Eq. ( |2.70[ ) and ( |2.72[ ) to mea- 
sure d. Indeed the two agree remarkably well - our pre- 
scription for the free energy of the dual theory passes this 
consistency test at zero temperature. 



3. Finite temperature 



The action Eq. (2.43) is only suitable at zero tempera- 



ture. We turn on a temperature in the little string theory 
by heating up the 5-branes. The resulting geometry is a 
black brane geometry 



g = L(-f dt * + dx 2 ) + — r 

■' R y J rRf 



rRg t 
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FIG. 4. The rescaled critical horizon radius as a function of 
a in the little string theory. Over this domain e 2 "^" decreases 
by a factor ~ 6.64 x 10 6 , while the rescaled radius is relatively 
constant. The critical radius therefore scales as c ~ e -27r//c \ 



where is the horizon radius and the dilaton and form 
fields are as before. This black brane has some interesting 
thermodynamics: its temperature is independent of r^, 



T = 



1 



2vri?' 



(3.43) 



This fact informs us that (i.) the little string theory has a 
Hagedorn phase and (ii.) the bulk theory at any (order 1) 
temperature is dual to that Hagedorn phase. Consistent 
with this result, the energy and entropy densities of the 
black brane grow with r^, 

(3.44) 



In both the bulk and field theories, entropy and energy 
can be added by increasing without raising the tem- 
perature. 

We will now measure how the BKT transition is mod- 
ified in the presence of this background entropy density. 
We will use a different parametrization for the embedding 
than in Eq. (2.32), writing the 3-sphere metric as 



3s3 



d9 2 



■6d<p 2 



] 6d^ 2 



(3.45) 



and letting 9 = 8(r) depend on the radial coordinate. 
As before our brane embeddings will be specified by <p = 
x 4 = x 5 =constant. The Legendre transformed brane 
action density now takes the form 



S 5 = -N J drv/l + r 2 /^' 2 ^ 



B 2 



r z cos^ 



(3.46) 

Before going on, we should pause to note an important re- 
sult. The linearized action for small fluctuations around 
the chirally symmetric 9 = is now 



£~ -Af\^p 2 + B 2 + r 2 



r 2 f6 n 



N{B 2 



l )B 2 



2\jp? + B 2 + r 2 

(3.47) 
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which, as for the ABJM theory, now exhibits three sep- 
arate limits at small horizon radius. At large r ^> fi, 6 
fluctuates like a stable scalar field and there is the AdS 2 
region at medium r/, <C r <C /i. However, there is now 
a near-horizon region r — <C 77, , conformal to Rindler 
space, where the field is stable. Following the same ap- 
proach as in Sec. |III A[ we describe this system with an 
equivalent quantum mechanics. This problem has a po- 
tential that has the form of an inverse square potential 
regulated at both small and large distances. As before 
a single tachyon will form as fj. is decreased, signalling a 
continuous transition rather than exponential scaling. 

Moreover we can estimate the scaling of the critical 
horizon radius at fixed a [23 . Using the same logic as in 
Sec. |III A| we find that to exponential accuracy the finite 
temperature transition occurs at 



n.. 



= -2tt/c 



(3.48) 
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FIG. 5. A log-log plot of the condensate as a function of 
magnetic field at nonzero temperature in the little string the- 
ory. The black dots are numerical data and the thin blue line 
is a fit with slope .500 ~ 1/2. This slope measures the critical 
exponent v, which takes the mean-field value. 



We confirm these results numerically by solving the 
embedding equations at nonzero temperature. We do 
so again by shooting inwards from the boundary. Our 
boundary condition there is that 9 is normalizable, for 
which we obtain a large— p series solution 

oo „ 
n=0 

where the 9 n are recursively determined by 9q. We use 
this series solution to set initial conditions for numerical 
integration. We then integrate down very close to the 
horizon, where we see if our solution satisfies the infrared 
boundary condition that 8 is regular at r^. 

We plotted some of our results in Sec. |T] We plotted 
the condensate in Fig. [lj unfortunately, we have suffered 
some systematic errors in our numerics for the free energy 
at nonzero temperature and so we do not plot them heraj 
We go on to present more data in this subsection. First, 
we plot the critical horizon radius as a function of a and 
find that it indeed scales exponentially in Fig. [3] Finally, 
we measure the static critical exponent 

\c\ ~ ( Mc (T) - At)", (3.50) 

for the transition at Th — 1CP 12 in Fig. [5] and find it to 
be mean-field, v = 0.500 ~ 1/2. 



4 Away from the chiral transition at nonzero temperature the em- 
beddings asymptote to the zero temperature result for r — 2> 
r/j. For small horizon radii the free energy should then also 
asymptote to the zero temperature data. We do not find this, 
implying some systematic errors in our determination of the free 
energy at finite temperature. 



IV. DISCUSSION 

A. IR AdS 2 

We have identified quantum holographic BKT transi- 
tions in both flavored ABJM and little string theories. 
The mechanism in the bulk that drives the transition is 
clear: in each case the field dual to the chiral symme- 
try breaking order parameter becomes unstable near the 
bottom of the flavor brane. The radial equation for the 
field in this region is that of a scalar in AdS 2 with a mass 
that depends on the control parameters. The holographic 
BKT transition occurs when the mass drops below the 
Breihtenlohner bound for stability in this effective 
AdS 2 region. 

We would like to use this information to infer proper- 
ties of the mechanism in the field theory that gives rise 
to the holographic BKT transition. As in we find it 
instructive to consider a multiplicity of systems that do 
not realize holographic BKT transitions. The flavored 
little string theory, in particular, is one of many systems 
whose holographic dual has supersymmetric flavor branes 
probing the near-horizon p-brane geometries [71] . These 
geometries are given by 

g = Z{fT? g P ,i + Z{rpr {dr 2 + r 2 g s s~ P ), 

2 

e*<xZ(r)^, Z(r)= T —, (4.1) 

where R is related to the coupling of the Yang-Mills the- 
ory on the p-branes and there are N units of (p + 2)- 
form flux. For p < 7 these field theories typically flow to 
flavored maximally supersymmetric Yang-Mills theories 
in limits that vary with p. Most supersymmetric probe 
branes in these backgrounds fall into one of three cate- 
gories: (i.) D(p + 4) branes that wrap an § 3 c S 8_p , dual 
to codimension-0 flavor, (ii.) D(p + 2) branes that wrap 
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an § 2 , dual to codimension-1 flavor, and (iii.) Dp branes 
that wrap an S , dual to codimension-2 flavor. The oth- 
ers have flavor in either (0 + 1) or (1 + 1) dimensions, for 
which wc cannot turn on a magnetic field. 

The probe brane action for the theory with 
condimcnsion-fc flavor at nonzero magnetic field and den- 
sity is not hard to obtain. Following the derivation of 
Eq. (2.381, we find that it is proportional to 



S P .k oc J cF +1 - k xdp^l^^d? + p2(3-fc) (i + J*L 

(4.2) 

where d is proportional to the constant of integration 
associated with A' in the DBI action, d£/dA' Q . For most 
cases this constant is interpreted as the baryon density of 
the dual theory. For others (like the flavored little string 
theory we consider), it is the chemical potential for the 
baryon density. Also, p + 1 — k must be greater than or 
equal to 3 so that we can add a magnetic field. 

The Lagrangian for embeddings close to the chirally 
symmetric one y = is 



C oc -yjd 2 + B 2 pP- 2k - 1 + p 2 (3- fe ) 

B 2 p-2k-3 



y_ 

2 

(7- P )y 2 



yjd 2 + B 2 pP- 2k ~ 1 +p 2 ( 3 - fe ) 4 



(4.3) 



This Lagrangian has an infrared limit that depends on 
whether p — 2k — 1 is positive, negative, or zero. They 
have the case structure 

1. p — 2k — 1 < 0: There is one such case, p = 4 and 
k = 2. It has a limit where y/p obeys the equation 
of motion of a m 2 = — 1 scalar in AdSa and so is 

2 

unstable - the symmetric embedding is unstable at 
all densities. 

2. p — 2k — 1 = 0: There are two cases, p = 3, k = 
1 and p = 5, k = 2. These are respectively the 
D3/D5 system and the flavored little string theory 
we consider in this work. In these cases y/p obeys 
the equation of motion in the infrared of a m 2 = 
— (3 — k)B 2 /(B 2 + d 2 ) scalar in AdS2, which passes 
below the BF bound for d < V7B in the D3/D5 
system and d < \/SB in the little string theory. In 
each case there is a holographic BKT transition at 
this magnetic field. 

3. p — 2k — 1 > 0: In all of these systems, y/p behaves 
in the infrared as massless scalar in AdS2- A large 
enough magnetic field will break chiral symmetry 
but the transition will not be of the holographic 
BKT type. From the numerical work done in the 
D3/D7 system [2T] we expect the transition to be 
second-order with mean-field exponents. The one 
exception is p = 6, k = where the magnetic field 
does not trigger chiral symmetry breaking at zero 
density; presumably there is no transition in that 
system at any magnetic field. 



The only theories in this class that realize holographic 
BKT transitions in the {B,d) plane are the D3/D5 and 
D5/D5 systems. We would like to understand what is 
special about these theories. In [22], we speculated that 
the special feature was that the theory lived in (2 + 1)- 
dimensions so that B and d had the same operator dimen- 
sion. The D5/D5 system is a striking counterexample, 
but can be understood similarly in terms of generalized 
conformal symmetry. 



B. Generalized conformal symmetry 

Most of these field theories run in the infrared (except- 
ing flavored J\f = 4 SYM). However, it is clear from their 
holographic duals that a generalized conformal symme- 
try emerges there [75] . Dilatation in the boundary the- 
ory is realized in the bulk by simultaneously scaling the 
field theory directions x and the radial coordinate. For a 
probe brane that adds codimension-fc flavor, the volume 
of the brane 



V k = J d^- 2k ^^P\g\ 



is invariant under 



X I— > A Q 



4(p + l-fc) 

r i— > A i-«+p(s-p) r . 



(4.4) 



(4.5) 



There is then a generalized dilatation invariance in the 
dual flavor, under which operators will generally assume 
a dimension that differs from its canonical one. As an 
example, let us obtain the dimensions of the magnetic 
field and chemical potential/density in these theories. 

To find these dimensions we need only consider the 
U(l) gauge field in the weak field limit, 



S EM = -N f T p+4 _ 2k J d» +5 - 2k £e-*^P[g} 



F2 



(4.6) 



This action is invariant under T>\ when the components 
of the field strength assume particular eigenvalues under 
T>\. For example, F 0r must have dimension and Fij 
dimension 2(p + 1 - fc)(7-p)/(l - 4fc +p(8 — p)). From 
this we learn the dimensions 



4(p + 1 - k) 



[d] = (3-k)[p] 



l-4k+p(8-p) : 

2(p+l-fc)(7-p) 



[B] 



l-4k+p(8-p) 



(4.7) 



where \p] comes from the scaling of the leading term of 
A and [d] from the near-boundary expansion of A a — 
p — d/J\f(2 — k)p 2 ~ k + .... These dimensions obviously 
differ from the canonical ones [p] — 1, [B] = 2. 

Now we note that the only (p, k) for which [d] = [B] 
are (3, 1) and (5, 2). These are the D3/D5 system and the 
flavored little string theory we study in this work. In each 
case B and d have dimension 2. Indeed, from our analysis 
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in Sec. IV A we know that these are the only systems in 
this class that realize holographic BKT transitions. 

What about the other (p, k) systems? It turns out 
that the one system with p — 2k — 1 < 0, (4,2), has 
[B] > [d] = The reader will recall that there is no 
chiral phase transition in this theory; chiral symmtrey is 
broken for all B ^ 0. Moreover, the other systems with 
p - 2k - 1 > have [B] < [d]. 

These results are extremely suggestive in terms of 
renormalization group flow. Consider a (generalized) 
conformal theory deformed with some set of control pa- 
rameters. These will induce flows of conjugate opera- 
tors. In particular, the infrared of the theory will be 
dominated by the flow of the lowest dimension operator 
which will be induced by the highest dimension control 
parameter. If there are two or more control parameters of 
the same dimension are turned on, they can dominate the 
infrared, provided that they remain the same dimension 
(relative to each other) after being turned on. In such 
a theory there will be non-trivial renormalization group 
flow as these operators compete to arbitrarily low energy 
scales. We claim that these conditions are sufficient to 
engineer a holographic BKT transition if one of these 
operators tends to restore a symmetry and the other to 
break it. At the least, this sort of scenario is markedly 
different than the usual Ginzburg-Landau- Wilson picture 
for phase transitions [86] . 

This picture is certainly consistent with all known 
quantum critical points embedded in string theory. The 
three probe brane systems with [B] — [d] realize holo- 
graphic BKT transitions. The theory with [B] > [d] = 
has an infrared dominated by the magnetic field and 
there is never a chiral transition. In the theories where 
[B] < [d], the infrared is density-dominated and the chi- 
ral transitions are presumably second-order. 

We have two loose ends to discuss. The first is the 
scalar mass operator in the flavored little string theory. 
How does it fit into this picture of generalized conformal 
invariance and BKT transitions? A simple computation 
shows that [M] = 2 in this theory, so that all three con- 
trol parameters p, B, and M have the same generalized 
conformal dimension. We are not surprised, then, to find 
that there is a line of BKT transitions in the (p, B, M) 
plane. Second, there are two more systems for which 
[fj] = [B]: (p,k) = (5,0) and (5,1). Why do wc not 
see holographic BKT transitions in these systems? Our 
answer is simple. In each case the dimension of the den- 
sity is larger than [pi\. As a result the deep infrared will 
be dominated by the density, rather than an interplay 
between the magnetic field and fi. 
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Appendix A: Computation of the spectrum 

We obtain the meson spectrum of the flavored ABJM 
theory by studying small fluctuations around the embed- 
ding y = m. There are four types of fluctuations corre- 
sponding to the four types of fields on the flavor brane, 
one corresponding to the transverse scalars y and those 
related by 50(3) x symmetry and three to the gauge field 
on the brane. We begin with a fluctuation of the trans- 
verse scalar y — m + Sy. 

The linearized Lagrangian for Sy is equivalent to 

duyd^y 



C 



5yl 



\J ' p 2 + m 2 



5y L 



(p 2 + m 2 ) 2 



(Al) 



where P[g] is the induced metric on the background y = 
m, di is a derivative with respect to a coordinate on C3 
and the i's are contracted with the metric on C 3 . By 
separability, the equation of motion for 5y can then be 
solved by an eigenfunction expansion. They are given by 



Sy k ,n,l = e l *- X <l> v> n{p)y 3 {Cs), 
where the y j (C3) are harmonics on C 3 , 

D i D i y l (c 3 ) = -\ j y j (c 3 ), 



(A2) 



(A3) 



and Di is the covariant derivative on C3. Since at any 
constant p C3 is a squashed MP 3 we can obtain the har- 
monics on C3 from the harmonics on MP 3 . These fall into 
(§, I) multiplets of the 50(3) x 50(3) isometry of MP 3 
for integer j/2. Under the pullback from the squashed 
RP 3 to the round one, these harmonics form a complete 
set on the squashed space. However, the reduced isome- 
try of the squashed MP 3 implies that these harmonics can 
mix and that their eigenvalues may be dependent upon 
the squashing parameter y 2 /p 2 . 

Actually, only the second 50(3) is broken to U(l) in 
the squashed space. The harmonics on C3 therefore fall 

into (l; |, ,j y^J 5 representations of 50(3) x [7(1) for inte- 
ger j/2 and m/2. There is only one such representation 
for each j/2, m'/2, so the scalar harmonics do not mix 
on C3. It is a short exercise to show that, however, the 

eigenvalues of the |, harmonics are modified to 
depend on y 2 /p 2 through 



j(j + 2) -m> [ 



(A4) 



5 This notation indicates both of the 50(3) quantum numbers as 
well as the U(l) charge. 
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The equation of motion for a mode Eq. ( A2 1 is then 

k 2 



P 



(p 2 + m 2 ) 



2\2 w 



1 

4p 



2 ^j,m'4>y 



= o, 

(A5) 



which has the solutions 



Kj,± = p ±! 



V 



,2x4/3 



-j ± ml + p 2 + j ± ml + (3 ml p 
4 ' 4 ' ± T' _ ^2 



(A6) 

2 



where we have defined j3 = l + yl — 4fc 2 /m 2 . We obtain 
the meson masses from these solutions by finding those 
—k 2 for which we satisfy boundary conditions at small 
and large p. Near the bottom of the brane at small p the 
4>y,± behave as 



(A7) 



Our small-p boundary condition is simply that 5y is reg- 

'j/,sgn(m')- At lar S e P 



ular so that the correct solution is 
where we have 

r(i±i)r(i±f) 



2+j±m'+/j 
4 



•0(p)" 



(A8) 



the proper boundary condition is that <j) y is normalizable. 
That is, the coefficient on the ppl 2 term must vanish. 
This can only be accomplished for particular k 2 so that 
the denominator is sitting on a pole. These are given by 



k; = — (2 + j+\m'\+in)(4 + j + \m'\+4n). (A9) 



These are the meson masses (-squared) of the dual theory. 
Moreover, the large- p behavior of the solutions Eq. (A8) 
yields the dimensions of the dual operators, A = 2 + |. 
These operators are of the form ^(AiBi)^ 2 ^, where the 
if> are the flavored fermions and the Aj's and B^s are 
bifundamental fields of the AB JM theory [65] • 



We move on to consider the fluctuations of the gauge 
fields. Following [S2HBS], these can be broken into three 
classes, 



Type I 
Type II 
Type III 



A„ 



A^ = 0, A p = e ik - x 4> n ,n(p)y j (C3), Ai=e ik - 
A^ = 0, A p = 0, A i = e tk - m <f> III , n (p)y i i (C 3 ), 



^//,„(p)A^'(C 3 ), 



(A10) 



where we have indexed the internal directions with i and 
the vector harmonics of spin j/2 under the unbroken 
SO(3) on C 3 are labeled by y? (C 3 ). These obey 



D z D*yi(c 3 ) - Kyf(c 3 ) = -A^(e 3 ), 
D'y? = o, 



(All) 



where [Di, Dk] = Rik is the Ricci curvature. 

We have been able to solve the meson spectrum for the 
Type I fluctuations. Due to the background Ramond- 
Ramond 3-form potential the gauge field obeys 



0. 



(A12) 



where we have chosen the orientation (tp,9,^>) to be pos- 
itive. For the Type I fluctuations the b = p, equations 
impose the usual condition £ • k = 0, the b = i equations 
are automatically satisfied, and the b = p equation is 



(gV ' p 2 + m 2 ^) 
y/ p 2 + m 2 



pk 2 



(P 2 



1 

Tp 



(A13) 



where the A,- m > are those in Eq. ( A4 1 . This is the same 



equation of motion as for the transverse scalar Sy and so 
has the same solutions and meson spectrum, 

9 

771 

-kj = — (2 + j + \m'\ + 4n)(4: + j + \m'\ + 4n). (A14) 

The dual operators have the dimensions 2 + | . 

Moving on to the Type II fluctuations, we find that 
the b = p equations relate 4>n & n d (\>n by 



(py/p 2 + m 2 4> II )' A ~ 



P 



\/ P 2 + 771 2 



4p 2 



5/7- 



(A15) 



Before going on, note that for j = this equation reduces 
to {py/ p 2 + mP^jx)' — 0, which has the solution 



9/7 oc 



p V p 2 + ■ 



(A16) 



which is singular at p = 0. We therefore consider only 
j > modes. Using the condition Eq. (A15), the b = p 
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and b — i equations are identical, 
I p^J p 2 + to 2 



4 m 2p2 £.2 

Aj, m '(™ 2 +p 2 ) 2 ) 



1 



J il 



(A17) 

where the Aj im < are as in Eq. (A4|. This equation of 
motion presumably has analytical solutions in terms of 
hypergeometric functions and an interesting fc-dependent 
prefactor. Unfortunately, we have not been able to find 
it. The best we have been able to do is to solve Eq. ( A17 ) 
at k = 0, 



_ i ' — j±m' 2 + j ± to' m 1 

v, =P ± ^ 2 F 1 [^- , J - ;1±- 

4 4 2 rn 



At large p these solutions go as 



/>// ~ + p 



2 

(A18) 



(A19) 



so that the dimensions of the dual operators are indeed 
2 + |. Nonetheless, the meson sp ectru m could be ob- 
tained by numerically solving Eq. (A17); we have done 



this for various j, to' and found a meson spectrum that 
agrees with the masses Eq. ( A9 ) to approximately one 
percent. 

The last set of modes, the Type III fluctuations, is 
by far the hardest to solve. To see this, recall that 
there are two types of vector harmonics on a round 
§ 3 and therefore on a round MP 3 « S 3 /Z 2 , those that 
fall into (^r 1 , ^f^ 1 ) representations and the gradients of 
the scalar harmonics, Diyi , which fall into (I, |) rep- 
resentations. There are then three representations that 



can mix for general j, to': 



(i 


j 




( j . J-2 


m' \ 




2 


2 ) 


^2' 2 


2 ) 



( j . j+2 m'\ 
\2 : 2 ' 2 J ' 



Since 



DiD i (D k y 3 ) = 



and 



(A20) 



there are only two linearly independent eigenmode rep- 
resentations for each j,m'. These modes will generally 
mix rather nastily. There are some exceptions. First, 



l. J+ 2 rnf+2 
2 ' 2 ' 2 



for j = there are only the ^0; 1, vectors, so these 
do not mix and are simple to analyze. Second, for the 
maximal \m'\ = j + 2 there is only the ( 
vector. Finally, the near-maximal case \m'\ = j yields a 
single mode which is a combination of the (|; I, |) and 

(§' ^1T' i) representations. 

It is easy to perform the harmonic analysis for the j = 
fields on the internal space. Unfortunately, we cannot 
solve the corresponding sectors of meson spectra. We 
look at these fields anyway in order to verify a claim in 



Sec. II A 2 about their near-boundary expansion. The 



eigenvalues are 

m! = ±2 : A = 4, m' = Q : A 
and the fields also obey 



p 2 



(A21) 



1, to' = ±2, 



d kl d k y? = -23>°Vfe < A, m' = 0. (A22) 



The equations of motion for these fields Eq. (A12) are 
then 



{p{p 2 +m 2 f/ 2 4>' IH )' (2p 2 -m 2 ) + ^ 



1 '1 - 



p(p 2 + m 2 ) 3 / 2 
for m' = ±2 and 



p 2 (p 2 + 



</>/// = o 



(p 2 +TO 2 ) 5 / 2 



(2p 2 + fc 2 )0 



77/ 



-\/p 2 + m : 

(A23) 

for to' = 0. In each case the large-p solution for <f>in is 
given by 



9/7/ 



ox 



«2 

n2 



(A24) 



As discussed in Sec. II A 2 the dimension of the dual 



operator in the supersymmetric theory is 1 so we should 
treat a 2 as a source. We do so in Sec. |IIC T] and find the 
one-point function of the dual operator. 
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